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SECTION  1 


INTRODUCTION  AND  SUMMARY 


A.  PROGRAM  OBJECTIVES 

The  primary  objectives  of  this  program  were  to  find  and  develop  new 
infrared  (IR)  transmitting  materials  and  provide  new  data  on  the  electro¬ 
optic  (EO)  properties  of  those  materials  most  likely  to  have  an  EO 
coefficient  an  order  of  magnitude  higher  than  those  of  materials  pres¬ 
ently  in  development  for  tunable  filters.  A  higher  EO  coefficient  would 
allow  an  engineering  margin  on  designs  for  increased  voltage  from  drive 
chips.  In  addition,  numerous  applications  can  be  envisioned  for  such 
materials  in  sensor  and  communication  systems.  Since  tensor  symmetry  is 
the  same  as  for  the  EO  tensor,  such  a  material  is  also  likely  to  have 
excellent  piezoelectric  properties  and  thereby  have  many  uses  as  an 
ultrasonic  transducer.  Many  of  the  materials  can  also  be  used  in 
devices  requiring  collinear  acoustooptic  effects. 

B.  BACKGROUND 

Approximately  three  years  ago,  personnel  at  Hughes  Research 
Laboratories  (HRL)  reviewed  the  requirements  for  an  IR  transmitting, 
spectrally  tunable  filter  for  the  HALO  program.1  We  concluded  that 
both  the  total  power  consumption  and  the  drive  power  density  would  be 
excessive  if  an  acoustooptic  filter  were  used.  This  led  us  toward 
the  invention  of  an  electronic  analog  to  the  A0  filter.  Such  a  filter 
would  require  orders  of  magnitude  less  drive  power.  Progress  on  the 
reduction  to  practice  of  this  invention  has  been  excellent,  but  develop¬ 
ment  work  is  not  without  technical  risks.  The  test-sample-acquisition 
and  materials-development  tasks  we  have  already  completed  in  pursuit 
of  the  filter  provided  us  with  a  data  base  of  EO  crystals  (mostly 
crystals  having  the  chalcopyrite  structure  (e.g.,  AgGa$2)»  which  was 
selected  for  the  HALO  program.  A  material  having  an  EO  coefficient  an 
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order  of  magnitude  greater  than  that  of  AsGaS2  would  reduce  required 

maximum  drive  voltages  to  less  than  100  V  maximum. 

2 

A  parallel  study  made  at  HRL  concerned  with  fundamental  EO  con¬ 
cepts  led  to  the  extension  of  the  number  of  allowable  crystal  classes 
suitable  for  application  to  tunable  filter  designs  when  crystallographic 
orientations  other  than  those  of  the  principal  axes  are  used  for  the 
selected  direction  of  light  propagation.  Figure  1  shows  the  crystal 
classes  as  indicated  by  the  point  groups,  which  are  useful  for  devices 
employing  the  longitudinal  EO  effect  using  principal  axes  (i.e.,  unro¬ 
tated)  and  those  additional  point  groups  that  become  applicable  when 
principal  axes  are  rotated.  Figure  2  shows  a  similar  extension  of  point 
groups  for  materials  having  the  transverse  EO  effect.  Materials  falling 
into  the  indicated  point  groups  are  applicable  to  tunable  filters 
recently  designed  at  HRL.  The  crystal  classes  containing  materials  use¬ 
ful  in  collinear  acoustooptic  devices  are  shown  in  Figure  3. 
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Figure  1.  Point  groups  for  longitudinal  EO  effect, 
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TYPE 

POINT  GROUP 

CUBIC 

23,  M3,  432,  43M.  M3M 

HEXAGONAL 

TRIGONAL 

(^T)  3.  [~32T|  (^)  3M 

TETRAGONAL 

(^4^)  C^)  4/M,  422.  (4MM^)  42M  4/M  MM 

- 

ORTHORHOMBIC 

222.  (^M2^  MMM 

~  MONOCLINIC 

Ci)  Ci) 2/M  _  — . 

TRICLINIC 

©  T 

TRANSVERSE  I - 1  TRANSVERSE 

UNROTATED  I - ^  ROTATED 


Figure  2.  Point  groups  for_ transverse  EO  effect 
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ACOUSTO  OPTIC  (COLLINEAR) 


Figure  3.  Point  groups  for  collinear  acoustooptic  effect. 

The  use  of  axis  rotation  considerably  widens  the  selection  of  new 
materials  available  for  tunable  filters.  However,  an  extensive  litera¬ 
ture  search  reveals  that  very  little  data  was  available  on  many  of  these 
materials. 


C.  SUMMARY 

We  have  succeeded  in  identifying  two  new  materials,  Cdln^Te^  and 
AgGaTe^,  which  have  the  potential  of  having  an  EO  constant  at  least  an 
order  of  magnitude  higher  than  those  of  materials  currently  available  for 
EO  tunable  filters. 

This  project  has  also  led  to  the  development  of  synthesis  tech¬ 
niques  for  materials  preparation  and  crystal  growth  of  binary  and  ternary 
chalcogenide  materials.  We  have  worked  out  a  technique  to  widen  the 
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range  of  materials  as  far  as  resistivity  is  concerned  for  measurement 
of  dielectric  constant  by  the  use  of  blocking  contacts. 

Our  predictions  for  new  materials  with  high  dielectric  constants 
and  hence  the  potential  for  high  EO  coefficients,  based  on  structural 
analogs,  was  proven  both  by  experiment  and  by  the  theoretical  work  of 
A.  Yariv  and  C.  Shih.  In  a  parallel  program,  Shih  and  Yariv  developed 
a  techhnique  for  calculating  EO  coefficient  based  on  the  localized  bond 
charge  model.  Their  results  on  binary  compounds  are  in  excellent  agree 
ment  with  experiment,  and  they  calculated  very  high  EO  coefficients 
for  the  materials  we  predicted. 

We  recommend  the  continuation  of  this  type  of  materials  study, 
specifically  to  obtain  single-crystal  samples  large  enough  to  make  EO 
coefficient  measurements  —  an  aspect  we  did  not  complete  during  this 
program. 
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SECTION  2 


MATERIALS  SELECTION 

This  investigation  involved  the  selection  of  new  materials  based  on 
compound  semiconductors  that  have  point  groups  appropriate  for  a  pro¬ 
grammable  spectral  filter.  The  prime  candidates  were  chosen  based  on  a 
semiconducting  analog  to  the  best  observed  EO  materials  in  the  visible 
spectrum  range  (i.e.  ,  LiNbO^  and  LiTaO^). 

We  applied  the  principles  of  crystal  chemistry  to  a  structural 
analog  of  lithium  niobate  within  an  allowed  crystal  class  to  find  a 
similarly  large  E-0  coefficient  in  a  chalcogenide  (IR  transmitting) 
material.  A  detailed  analysis  of  the  lithium  niobate  structure 
(Figure  4)  shows  that  it  consists  of  a  planar  sheet,  structure  of  anions 
(oxygen)  in  approximately  hexagonal  close  packing,  which  results  in 
distorted,  partially  occupied  octahedral  sites  —  essentially  a  defect 
structure.  These  sites  are  one-third  occupied  by  niobium  and  one-third 
by  lithium;  the  remainder  are  vacant.  The  sequence  of  distorted 
octahedra  along  the  c-axis  with  Nb  at  the  origin  is  Nb-V-Li-Nb-V-Li, 
where  V  is  a  vacancy.’ 

A  plot  of  the  vertical  position  of  Nb  (in  an  octahedral  site) 
versus  potential  energy  (Figure  5)  shows  curves  for  temperatures  above 
and  below  the  Curie  point  (T^).  The  existence  of  curve  2  (with 

T  -  T  )  indicates  a  large  ionic  contribution  to  the  polarizability, 

c 

contributing  to  a  large  EO  coefficient,  which,  for  LiNbO  ,  is  32  x 

-12  J 
10  m/V. 

An  obstacle  in  seeking  a  semiconductor  analog  appeared  to  be  the 
tetrahedral  coordination  typical  of  semiconducting  materials  (e.g., 
zincblende  and  chalcopyrite  structures).  However,  our  study  uncovered 
a  group  of  materials  of  the  general  formula  that  are 

closely  analogous  to  LiNbO  and  thereby  show  considerable  promise  of 

_  3 

having  high  EO  coefficients.  These  materials  have  been  described  as 
crystallizing  in  the  crystal  classes  3m  (same  as  LiNbO^),  4,  and  42m 
(same  as  chalcopyrites) ;  all  of  these  are  potentially  usable  for  tunable 
filters. 


The  3m  point  group,  of  which  ZnIn2S/  is  a  member  (Figure  6),  has  a 
unit  cell  consisting  of  12  closely  packed  sulfur  layers.  Octahedral 
sites  contain  one-half  the  total  number  of  In  atoms;  Zn  atoms  and  the 
other  half  of  the  In  atoms  are  in  tetrahedral  sites.  An  inverse  struc¬ 
ture  has  also  been  described  in  which  the  divalent  and  one-half  of  the 
trivalent  atoms  are  in  octahedral  sites;  one-half  of  the  trivalent  are 
in  tetrahedral  sites,  leaving  a  fraction  of  vacancies  in  each  type  of 
site.  The  inverse  structure  has  been  described  for  Xln^S^,  where  X  is 
Mg,  Fe,  Co,  Ni,  CdGa^S^,  or  Cdln^Te^.  A  structural  comparison  of  two 
chalcogenides  with  LiNbO^  is  shown  in  Table  1.  Existing  theories  of 
the  EO  effect  in  crystals  predict  that  the  EO  coefficients  of  these 
materials  should  be  significantly  higher  than  those  of  the  chalcopyrites . 
There  is  little  data  in  the  literature  concerning  the  preparation  and 
properties  of  these  materials. 

Because  the  lead  time  required  to  synthesize  single-crystal  samples 
is  long,  our  experimental  approach  involved  synthesizing  small  poly¬ 
crystalline  samples  of  a  wide  selection  of  candidate  materials  and  then 
measuring  the  dielectric  constant  of  each  at  ambient  and  liquid  nitrogen 
temperatures.  The  low-frequency  dielectric  constant  is  essentially  a 
measure  of  lattice  polarizability.  A  large  lattice  contribution  to  the 
EO  coefficient  is  essential  but  not  sufficient  for  a  large  EO  effect. 
Thus,  measuring  the  dielectric  constant  is  an  excellent  screening  method 
in  a  materials  survey  of  this  nature. 

A  second  phase  planned  for  this  program  was  aimed  at  (1)  growing 
single  crystals  of  at  least  two  materials  such  as  Znln^S^  and  Cdln^Te^ 
or  the  best  (highest  dielectric  constant)  from  the  materials  surveyed 
in  the  first  phase  of  investigation  and  (2)  subsequently  evaluating 
the  electrical  and  optical  properties  of  these  materials,  including  EO 
coefficient,  transmission  range,  resistivity,  refractive  index,  and 
birefringence. 
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Table  1.  Structural  Comparisons  of  Chalcogenides  with  LiNbO^ 
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Figure  6.  Structure  of  ZnIn„S,. 
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SECTION  3 


MATERIALS  SYNTHESIS 

A,  BINARY  COMPOUNDS 
1.  Synthesis  of  GeS9 

Our  first  attempt  to  synthesize  GeS_  (orthorhombic,  MM2),  we  used 

^  A 

a  reaction  procedure  similar  to  one  developed  at  HRL  for  AgGaS0.  High- 
purity  germanium  (greater  than  99.9999%  pure)  was  placed  in  a  fused 
quartz  boat  and  subsequently  in  a  sealed,  evacuated  fused  silica  tube 
with  the  required  amount  of  sulfur  (minimum  99.9995+%  pure)  to  react 
according  to  the  relation: 

Ge(s)  +  S_  (v)  GeS  (s) 

1  2 

at  temperatures  where  sulfur  vapor  (v)  reacts  with  solid  (s)  germanium 
to  form  the  solid  compound  GeS2»  A  schematic  of  the  apparatus  used  is 
shown  in  Figure  7. 

The  reaction  takes  place  readily  as  the  tube  is  moved  into  the  hot 
furnace  even  though  Ge  is  solid  (melting  point  =  956 °C).  However,  this 
procedure  yielded  two  products:  GeS0  and  GeS.  Although  the  experiment 
was  designed  to  prevent  formation  of  GeS  by  using  a  temperature  at 
which  that  compound  is  unstable  and  decomposes,  a  signficant  amount  of 
it  transported  to  the  cold  portion  of  the  tube  and  remained  intact. 

GeS2,  which  also  formed,  sublimed  as  well  to  the  cooler  portion  of  the 
tube.  Before  completion  of  the  run,  the  combined  sublimate  plugged  the 
tube,  isolating  the  remaining  elements  and  preventing  further  reaction. 
At  this  point,  the  run  was  aborted. 

To  avoid  a  repeat  of  this  problem,  we  modified  our  approach  using  a 
zone  refiner  (Figure  8)  with  two  furnaces  (two  zones).  In  addition,  a 
heated  outer  tube  employing  a  flat  nichrome  ribbon  to  maintain  a  mini¬ 
mum  temperature  throughout  the  tube  of  '\'400°C  was  established  to 
maintain  a  constant  overpressure  of  sulfur  above  the  germanium.  The 
zone  furnaces  were  heated  to  850°C  and  800°C,  respectively;  the  tempera¬ 
ture  of  these  furnaces  is  not  critical.  Rapid  reaction  was  observed 
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FURNACE 


T  -  850°C 


Figure  8.  Zone  refiner  used  for  formation  of  GeS^- 
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between  the  sulfur  vapor  and  germanium,  forming  GeS^  (white  in  color), 
which  sublimed  to  a  cooler  portion  of  the  tube.  A  total  of  six  zone 
passes  (three  passes  times  two  furnaces)  was  made  following  in  situ 
heating  during  which  the  germanium  sat  in  the  higher  temperature  furnace 
until  the  reaction  with  sulfur  appeared  to  have  gone  to  completion.  The 
resulting  product  was  white  crystalline  GeS^  with  a  small  residue  of 
yellow  amorphous  material  (GeS^).  The  crystalline  material  was  identi¬ 
fied  by  X-ray  diffraction  analysis  to  be  orthorhombic  GeS„  (see 

5  1 

Figure  9).  The  melting  point  of  GeS_  is  reported  to  be  840  ±  5°C. 

5  1 

The  equilibrium  phase  diagram  for  the  system  Ge-S  is  shown  in  Figure  10. 

Two  approaches  were  taken  to  grow  single  crystals  of  GeS^i  subli¬ 
mation  and  Bridgman  techniques.  In  the  preparation  of  GeS^  by  direct 
reaction  of  the  elements,  the  observation  that  GeS^  sublimed  readily 
presented  the  possibility  that  GeS2  could  be  grown  by  sublimation.  The 
compound  was  encapsulated  under  vacuum  into  a  1-cm  (i.d.)  tube,  20  cm  in 
length,  and  placed  in  a  furnace  with  a  flat  profile  at  700°C.  The  tube 
was  pulled  through  the  profile  at  a  rate  of  0.1  cm/hr  for  15  days.  No 
sublimation  took  place.  The  material  was  solid  and  polycrystalline. 

No  information  could  be  found  on  the  vapor  pressure  of  Ge$2  at  any 
temperature.  Since  the  compound  was  observed  to  sublime  and  possibly 
to  decompose  at  high  temperatures,  attempts  at  Bridgeman  growth  would 
increase  the  possibility  of  an  explosion  due  to  excessively  high  vapor 
pressures.  The  Bridgman  method  was  tried  in  a  furnace  set  at  850°C. 

The  rate  of  growth  was  2  mm/hr.  The  run  was  taken  out  of  the  furnace 
hot  before  the  last  remaining  liquid  had  frozen.  A  dark  brown  vapor  was 
observed  above  the  compound,  indicating  a  sulfur  pressure  over  the  melt 
and  hence  some  decomposition.  The  result  was  a  yellowish  crystalline 
material  identifed  as  Ge$2  by  X-ray  diffraction  (Figure  11),  from  which 
single-crystal  layers  could  be  peeled.  An  amber-colored  amorphous  sec¬ 
tion  (the  last  to  freeze)  was  also  identified  as  GeS2  by  microprobe 
analysis. 
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Figure  9. 
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X-ray  diffraction  pattern:  orthorhombic  GeS 
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Figure  11.  X-ray  diffraction  pattern  (indexed)  of  GeS2* 
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2.  Preparation  of  GeSe2 


GeSe^  was  synthesized  by  mixing  stoichiometric  amounts  of  Ge  and  Se 
in  two  vitreous  silica  ampoules  under  a  flow  of  nitrogen.  An  excess 
of  1.0  g  of  Se  was  added  to  provide  1  atm  of  Se^  vapor  at  the  processing 
temperature.  The  tube  was  then  evacuated  and  sealed. 

The  setup  was  heated  to  410°C  over  a  period  of  16  hr,  and  690°C  for 
30  hr  (boiling  point  of  Se  is  688°C).  The  color  of  the  Se^  vapor  was 
used  as  an  indicator  of  the  extent  of  completion  of  the  reaction.  When 
the  color  of  the  vapor  appeared  to  become  lighter,  the  temperature  was 
raised  to  725°C  for  18  hr,  then  to  740°C  (melting  point  of  GeSe^)  for 
four  days.  After  this  time,  no  Se^  vapor  was  observed  and  the  material 
sublimed  to  the  colder  ends  of  the  tube.  Part  of  the  tube  was  pushed 
out  of  the  furnaces  to  attempt  to  nucleate  and  grow  crystals.  The 
temperature  of  the  furnace  was  increased  to  765°C  in  two  days,  after 
which  the  setup  was  slowly  cooled  to  room  temperature. 

The  material  obtained  was  a  mica-like  mixture  of  gray,  red,  and 
yellow.  However,  X-ray  analysis  showed  the  material  to  be  GeSe2»  The 
color  differences  were  caused  by  variations  in  thickness. 

The  material  was  transferred  into  a  vitrious  SiO^  crystal  growth 
tube,  evacuated,  and  sealed.  Crystal  growth  by  the  vertical  Bridgman 
method  was  attempted  at  800°C.  The  resulting  materials  was  polycrystal¬ 
line.  Samples  were  sliced  and  submitted  for  dielectric  constant 
measurements. 

3.  Preparation  of  NbGe2 

The  preparation  of  NbGe^  was  done  using  the  constituent  elements 
powdered  into  an  intimate  mixture.  The  mixture  was  placed  in  a 
vitreous  carbon  boat  inside  an  aluminum  tube.  Under  a  nitrogen  atmos¬ 
phere,  the  setup  was  heated  to  1400°C  for  1.5  hr  and  1450°  for  2  hr. 

The  material  was  heated  to  1400°C  for  1.5  hr  and  1450°  for  2  hr.  The 
material  was  sintered.  X-ray  analysis  showed  the  material  to  be  a 
mixture  of  NbGe^  and  Ge. 


23 


B.  TERNARY  CHALCOPYRITES 

1.  AgGaTe2 

Polycrystalline  ingots  of  AgGaTe2  were  grown  from  the  melt  using 
the  Bridgman-Stockbarger  technique.  AgCaTe^  Is  believed  to  be  a  con- 
gruently  melting  chalcopyrite,  but  the  existence  region  has  not  been 
studied  in  detail.  Successful  synthesis  of  AgGaTe^  was  achieved  from 
Ga2Te^-rich  solution.  Our  measurements,  as  discussed  below,  indicate 
the  possibility  of  variations  in  stoichiometry  within  the  ingot. 

2.  Synthesis  of  LiGaS2 

Because  L^S  is  deliquescent  and  highly  reactive  to  water  in  the 
atmosphere,  the  ternary  chalcogenide  LiGaS^  was  prepared  by  reaction  of 
with  LiGaO^  rather  than  with  the  binary  compounds  Li^S  and  Ga^S^. 
Gallium  metal  was  first  dissolved  in  15  N  HNO^;  after  complete  dissolu¬ 
tion,  a  stoichiomteric  amount  of  LiNO^  was  added.  The  mixture  was 
evaporated  to  dryness  on  a  hot  plate  to  obtain  LiGa(NO^)^.  This  mater¬ 
ial  was  transferred  into  a  fused  silica  boat  in  a  fused  silica  ampoule. 

Since  the  physical  and  chemical  properties  of  LiGaS^  are  not 
reported  in  the  handbooks,  the  preparation  had  to  be  carried  out  cau¬ 
tiously;  a  great  part  of  it  depended  on  observation.  The  set-up  was  - 
heated  slowly  to  approximately  600°C  under  N2  to  decompose  the  nitrate 
to  file  oxide.  When  brown  fumes  of  nitrogen  oxide  were  no  longer  visible, 
the  temperature  was  raised  to  800°C  over  a  period  of  4  hr  under  l^S. 
Emerging  H^S  gas  flowed  into  an  outlet  connected  to  two  bubblers  contain¬ 
ing  NH^OH;  these  acted  as  traps.  X-ray  diffraction  showed  that  the 

material  was  mostly  LiGaO^. 

The  LiGaO^  was  ground  and  reheated  to  1000°C  for  2  hr.  This  time, 
te  material  was  heterogeneous:  colored  pink  and  white.  X-ray  analysis 
of  both  phases  gave  the  same  lines,  but  the  pink  phase  had  better  defined 
lines  and  the  LiGaO^  lines  were  not  as  intense.  There  were  other  lines 
present  and  these  could  not  be  identified.  (There  is  no  ASTM  index  card 
for  LiGaS^.)  These  lines  could  not  be  attributed  to  Li^S,  Ga^S^,  LiGa 
silicates,  Li  silicates,  or  Ga  oxides. 


24 


The  material  was  ground  again  and  then  reheated  to  1000°C  for 
3  hr  in  H^S  and  N^.  Melting  was  observed.  X-ray  analysis  showed  the 
same  lines  as  seen  during  the  previous  run.  Since  there  was  some  attack 
on  the  SiC^  (probably  due  to  the  Li^O),  the  run  was  repeated  using  a 
graphite  boat  at  1000°C  for  3  hr  under  K^S  and  N^.  The  final  product, 
which  was  mostly  LiGaS^,  still  retained  a  significant  amount  of  oxide. 

C.  DEFECT  CHALCOPYRITE ' COMPOUNDS 

TI  III  IV 

Compounds  of  the  type  A  are  known  as  defect  chalcopyrite 

compounds.  The  synthesis  of  these  compounds  (e.g.,  Znln^S^)  was  ini¬ 
tially  approached  by  vapor  transport  (using  iodine  as  a  transport  medium) 

because  of  the  various  properties  of  these  compounds,  including  high 

3 

melting  points  and  incogruent  melting  properties.  The  selection  was 
based  on  reports  of  crystal  growth  in  the  Russian  literature  as  cited 
in  Ref.  3.  We  used  mixtures  of  II-VI  and  III-VI  compounds  as  starting 
materials  for  the  reaction.  High-purity  ZnS,  ZnSe,  and  CdS  were  pur¬ 
chased  from  General  Electric  Company  and  further  purified  by  vacuum  sub¬ 
limation.  In  this  process,  a  dynamic  vacuum  is  used  (Figure  12);  pump¬ 
ing  is  continuous  through  a  liquid  nitrogen  trap  to  prevent  back- 
diffusion  as  well  as  to  enhance  the  removal  of  impurities  by  the  addi¬ 
tion  of  a  cryogenic  pump.  The  fused  quartz  tube  is  transported  (raised) 
through  a  thermal  gradient  as  indicated  in  Figure  12.  The  material  sub¬ 
limes  to  a  cooler  portion  of  the  tube;  many  more  volatile  impurities 
are  removed  and  trapped  by  pumping  while  others  segregate  according  to 
their  vapor-solid  equilibrium  temperatures. 

The  synthesis  of  inv°lved  a  two-temperature  process  using 

two  vitreous  silica  boats  in  a  vitreous  silica  tube  (Figure  13).  A 
stoichiometric  equivalent  of  sulfur  plus  a  small  excess  to  provide  1  atm 
overpressure  of  at  the  reaction  temperature  was  Introduced  into  the 
first  boat  at  one  end  of  the  tube  and  then  melted,  cooled,  and  pumped  to 
remove  occluded  H^O.  The  stoichiometric  weight  of  indium  was  placed 
into  the  second  boat.  Since  the  In  pellets  appeared  slightly  oxidized, 
the  In  end  was  heated  to  1000°C  for  2  hr  under  a  flow  of  keeping  the 
other  (sulfur-containing)  end  cool.  The  tube  was  cooled,  evacuated  to 
approximately  100  pm,  and  then  sealed. 
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Figure  12.  Preparation  of  In^Se^  or  In^S^, 
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Figure  13.  Sublimation  purification  of  II-VI  compounds. 

Precaution  was  taken  to  avoid  a  runaway  exothermic  reaction.  The 
preparation  involved  slow  heating  to  900°C  for  In  and  450°C  for  S  within 
a  12-hr  period.  The  last  step  consisted  of  heating  to  transport  all 
unreacted  sulfur  to  the  sulfur-containing  end.  The  bulk  material  pro¬ 
duced  was  identified  by  X-ray  diffraction  analysis  to  be  6  -  In^S^. 

Mixtures  of  binary  chalcogenides  can  be  used  as  source  material 
in  the  gas  transport  reaction  method  to  synthesize  ternary  compound 
crystals.  Following  the  preparation  of  ZnS,  CdS,  CdSe,  CdTe,  In^Se^, 
Ga2S^,  In^S^,  Ga2Se^,  and  In2Te^,  we  proceeded  to  investigate  the 
synthesis  of  the  following  tenary  compounds  based  on  the  reactions: 
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CdTe  +  In0Te,  +  Cdln.Te.  . 

2  3  2  4.- 

-  -  / 

1.  Preparation  of  Ternary  Compounds  by  Chemical  Vapor  Transport 

Chemical  vapor  transport  (CVT)  is  a  useful  technique  for  materials 
that  have  a  high  melting  point,  dissociate  at  the  melting  point,  or  melt 
under  elevated  pressures.  By  using  CVT,  lower  temperatures  can  be  used 
since  the  volatile  intermediate  phases  have  higher  vapor  pressures 
relative  to  the  material  of  interest.  In  addition,  this  method  allows 
separation  of  phases,  since  the  solid  condenses  within  the  section  of 
the  container  that  is  at  the  proper  temperature.  Compounds,  that  melt 
incongruently  (i.e.,  with  decomposition)  can  also  be  synthesized  since 
the  growth  takes  place  at  temperatures  below  the  melting  point  of  the 
compound.  Many  of  the  defect  chalcopyrites  of  the  type  fall 

into  this  category. 

The  mass  transport  rate  for  this  process  is  a  function  of: 

C>  Tube  geometry.  Increasing  the  cross  section  of  the 
tube  increases  the  rate  (which  is  controlled  by 
thermal  convective  processes);  increasing  the  length 
of  the  tube  slows  down  the  rate  of  transport. 
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The  concentration  of  the  transport.  At  low  concen¬ 
trations,  the  rate  is  slow.  As  the  concentration 
is  increased,  the  rate  increases,  slowly  at  first 
and  then  more  rapidly. 


The  average  temperature  and  the  temperature 
difference  (AT)  between  the  source  material  and 
the  growing  crystal.  The  temperature  difference 
will  be  of  minor  importance  whan  thermal  convective 
processes  are  dominant.  Higher  temperatures 
increase  both  vapor  concentration  and  convection, 
as  discussed  above.  A  greater  thermal  differential 
increases  the  transport  rate  since  vapor  pressure 
will  tend  to  equilibrate  to  the  temperature  at  the 
source,  while  condensation  occurs  in  the  cooler 


Iodine  as  Transport  Medium 


Initial  results  in  vapor  transport  reaction,  using  II-VI  and 

III-VI  binary  compounds  as  starting  materials  followed  by  the  addition 

of  iodine  into  the  ampoule,  were  not  universally  successful  (as  reported 
3 

in  the  literature  )  in  forming  the  desired  ternary  compound  (i.e., 

AIInIII„IV,  .  ,  .  .  ,  ,  '  . 

A  B2  )  or  even  in  completely  transporting  material  to  form  an  ingot 

(as  previously  described).  Those  compounds  that  had  not  previously 
yielded  usable  samples  for  evaluation  (either  by  reaction  and  recrystal¬ 
lization  in  situ  or  by  vapor  transport)  were  repeated  using  Cd^  and  Znl^ 
as  starting  materials.  In  these  cases,  iodine  is  still  the  transport 
medium,  but  advantages  are  gained  by  (1)  minimizing  the  loss  of  the  very 
volatile  elemental  iodine  in  the  preparation  stages  and  (2)  having 
constituent  elements  in  a  more  volatile  form,  e.g. ,  the  iodide,  which 
yields  a  highly  reactive  metal  atom  on  decomposition. 


(1)  ZnIn2S^  and  CdIn2S^  —  The  CVT  method  using  iodine  as  the 
transport  medium  initially  was  used  to  grow  Znln^S.  and  CdIn„S. ;  this 
technique  has  been  described  in  the  literature.  Two  tube  sizes  were 
used  for  growing  these  compounds:  each  had  a  1-cm  i.d.,  but  one  was 
20  cm  long  and  the  other  was  33  cm  long.  The  20-cm  tubes  were  pulled 
through  a  sloping  temperature  gradient,  while  the  33-cm  tubes  were  placed 
in  a  thermal  gradient  and  remained  stationary.  A  sloping  gradient  is 
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useful  along  the  crystallization  tube  to  find  the  optimal  temperature  for 
crystal  deposition.  The  tubes,  which  contained  ^25  g  of  crystal  source 
material,  were  loaded  with  iodine  in  various  amounts  and  sealed  off 
under  vacuum.  To  prevent  the  iodine  from  subliming  during  evacuation, 
the  tube  was  immersed  in  liquid  nitrogen. 

Five  runs  were  made  to  grow  Znln^S^.  Three  of  the  runs  consisted 

of  pulling  the  tube  through  the  temperature  gradient  from  930°C  to  700°C 

at  a  speed  of  2  mm/hr.  In  the  first  of  these,  the  tube  exploded;  in 

the  remaining  two,  the  runs  were  completed  without  incident.  The  other 

two  runs  of  the  five  were  done  in  situ.  In  those  two  runs,  the  tubes 

were  loaded  as  before  with  equal  molar  amounts  of  ZnS  and  In^S^,  iodine 

was  introduced,  and  the  tube  was  evacuated  and  sealed.  These  tubes 

were  placed  in  a  temperature  gradient  of  1000°C  to  700°C.  The  growth 

time  for  each  was  ''-'150  hr.  The  concentration  of  iodine  ranged  from  0.6 
3 

to  6  mg/cm  of  tube  volume. 

In  all  cases  described  above,  there  was  insignificant  transport. 

We  identified  Zn^In^S^  by  X-ray  analysis  (Figure  14)  as  the  resulting 
compound  that  dominated  in  all  the  runs.  A  small  amount  of  transported 
material  was  identified  as  being  Znln^S^  (Figure  15).  The  bulk  of  the 
material  sintered  and  underwent  grain  growth;  this  was  identified  as  being 
primarily  Zn^In^S^. 

Equal  molar  amounts  of  CdS  and  In  S_  along  with  iodine  equivalent 
3  ^  J 

to  6  mg/cm  of  tube  volume  were  sealed  under  vacuum  into  a  33-cm  tube. 
This  tube  was  placed  in  a  temperature  gradient  of  1000°C  to  700°C. 

After  150  hr,  no  significant  transport  had  taken  place.  -  The  material 
sintered  and  underwent  grain  growth.  X-ray  diffraction  analysis  was 
used  to  identify  it  as  CdIn2S^  (Figure  16).  A  1  cm  x  2  mm  sample  has 
been  cut  for  measurements. 

The  reason  for  the  lack  of  transport  in  the  above  experiments  is 
not  obvious.  The  tubes  were  of  sufficient  diameter  and  not  too  long; 
the  sizes  are  comparable  to  the  sizes  reported  in  the  literature.  Also, 
tubes  of  this  size  had  been  used  on  numerous  occasions  in  this  labora¬ 
tory  to  grow  binary  compounds  by  CVT  techniques.  The  concentration  of 
iodine  is* sufficient  without  getting  into  the  dangerous  situation  of 
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Figure  14 


X-ray  diffraction  pattern  of  Zn.In„S, 
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Figure  16.  X-ray  diffraction  pattern  of  Cdln^S^. 
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pressure  buildup.  The  temperature  difference  (AT)  is  large  enough  for 
transport  even  though  it  may  not  be  operative  since  convective  processes 
may  be  dominant.  The  length  of  time  in  the  furnace  was  certainly  ade¬ 
quate  for  transfer  to  occur. 

This  procedure  was  followed,  to  synthesize  Cdln^S^,  Cdl^Se^,  and 
ZnIn2S^,  using  1-cm-diameter  ampoules  and  calculating  iodine  concentra¬ 
tion  to  obtain  1  atm  of  iodine  at  the  run  temperature;  stoichiometry  was 
balanced  using  appropriate  II-VI  binary  compounds.  Runs  typically  were 
carried  out  around  1000°C  at  the  source  with  a  large  thermal  gradient  on 
the  order  of  400°C. 

Successful  results  in  transporting  and  forming  CdIn2Se^  were 
obtained;  the  compound  was  identified  by  X-ray  analysis  (Figure  17).  A 
small  amount  of  ZnS  +  In9S^  transported,  but  not  enough  for  identifi¬ 
cation.  No  transport  of  constituents  was  observed  in  the  CdIn2S^ 
experiment. 


50  40  30 
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Figure  17.  X-ray  diffraction  pattern  of  Cdln^Se.  , 

l  4 
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HC1  as  Transport  Medium 


Early  observations  of  little  or  no  transport  using  iodine  as 

a  carrier  gas  led  to  the  search  for  a  more  efficient  carrier.  Based  on 

optimistic  reports  in  the  literature,  we  selected  HC1  for  this  carrier 

gas.  Special  equipment  for  loading  this  corrosive  gas  into  the  ampoule 

had  to  be  acquired.  The  only  additional  significant  change  in  these 

runs  was  to  increase  the  tube  diameter  from  1  cm  to  2  cm  in  accordance 

with  theories  that  the  transport  rate  increases  exponentially  with 

increased  ampoule  diameter  (as  discussed  above).  Using  HC1  as  the 

transport  medium  at  a  pressure  calculated  to  be  1  atm  at  run  temperature, 

CdIn2S^  was  prepared  and  confirmed  by  X-ray  diffraction  analysis. 

Initial  runs  using  HC1  for  the  preparation  of  Znln^S^  indicated 

that  there  was  a  substantial  amount  of  the  desired  ternary  compound 

present  as  well  as  ZnnIn0S01  and  Zn„In0S,.  A  ternary  compound  close 

yoZJL  j/o 

to  CdGa^S^,  as  identified  by  X-ray  analysis,  was  also  formed  by  HC1 
vapor  transport  reactions.  Some  free  Ga^S^  (a  volatile  component)  was 
collected  from  the  run  and  identified.  In  addition,  the  formation  of 
CdGa^Se^  from  its  binary  constituents  (CdSe  and  Ga^Se^)  was  accomplished 

using  this  procedure. . -  - - - - - — - 

The  results  of  all  vapor  transport  reaction  studies  are  summarized 
in  Table  2. 

(1)  ZnGa2S^  and  CdGa2Se(!(  Synthesis  Using  HC1  as  Transport 
Medium  —  In  the  case  of  ZnGa^S^,  all  starting  material 
was  transported  to  the  conical  end  of  the  sealed  ampoule  at  tempera- 
tures  as  high  as  1100°C.  The  ampoule  contained  three  distinct  bands:  _ 
yellow,  white,  and  translucent.  Initial  DTA  investigation  of  ZnGa^S^ 
indicated  a  melting  point  in  excess  of  2100°C  (limit  of  current  DTA 
apparatus)  and  a  possible  phase  transition  above  1000°C.  Attempts  at 
single-crystal  growth  by  vapor  were  not  successful.  — 

The  sample  of  CdGa^Se^  showed  almost  no  transport,  yielding  only  _ 
a  thin  coating  showing  two  different  colors.  The  temperature  in  the 
hot  zone  was  900°C;  1  atm  of  HC1  was  maintained  at  the  operating 
temperature. 

Results  indicated  the  need  for  an  increased  temperature  (1000°C 
appears  useful)  and  for  an  increased  pressure  of  HC1  at  operating 
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Table  2.  Vapor 


Desired 

Compound 

Starting 

Materials 

Temperature 
Range,  "C 

Carrier 

Tube 

Diameter, 

cm 

Cd In  Sc. 

2  U 

CdSe  +  In2Se^ 

990  to  550 

*2 

1 

Cd  In  S. 

2  U 

CdS  +  ln,S3 

990  to  550 

l2 

1 

2  U 

ZnS  +  In^S^ 

990  to  550 

l2 

1 

7.nln,S4 

ZnS  -4- 

1010  to  625 

HC1 

2 

ttJIn-S. 

2  U 

CdS  +  In2S3 

1010  to  625 

HC1 

2 

CdS  +  Ca2S3 

650 

HC1 

2 

CdC.t-Sc* 

2  U 

CdSe  +  Ca^Se^ 

900 

HC1 

2 

Reaction  Studies 


Amount  of 
Source  Left 


Transport 


Remarks 


None 

All 

Much 


Much 


Yes 

Yes 

No 

— 

Little 

Two  sections  of  tube  coated: 
yellow  crystalline  In  high 
temperature  section;  red, 
yellow,  white  coatings  in  low 
temperature  region. 

Yes 

Yes 

(Black) 

Also  Zn^IngS2l  and  ZnjIn^S^ 

i 

Yes 

(Black) 

Yes 


Yes 


temperatures.  A  review  of  recent  crystral  growth  literature  uncovered 
a  report  of  the  use  of  2  to  3  atm  pressure  of  HC1  under  operating  condi¬ 
tions  for  a  successful  yield  of  ternary  chalcogenide  crystals  by  vapor 
transport. 

(2)  Preparation  of  ZnIn2S^  and  CdIn2S^  —  For  the  preparation 
of  Znln^S^,  the  starting  materials  were  the  respective  compounds  encapsu¬ 
lated  in  a  silica  tube  with  380  mm  HC1  at  ^25°C.  To  clean  its  tip,  the 
tube  was  placed  vertically  in  a  hot  furnace  for  one  day  so  that  the 
collection  end  of  the  tube  was  hotter  than  the  source  end.  The  tempera¬ 
ture  was  then  reversed,  and  the  tube  was  pulled  slowly  (0.8  mm/hr) 
through  the  temperature  profile  until  it  was  several  centimters  out  of 
the  furnace.  For  the  zinc  compound,  the  maximum  temperature  was  1000°C 
and  the  "growth"  time  was  18  days;  for  the  cadmium  compound,  the  maxi¬ 
mum  temperature  was  960°C  and  the  "growth"  time  was  19  days.  After  the 
run,  the  tube  was  air  quenched  to  room  temperature. 

In  previous  runs  with  HC1,  the  pressure  of  HC1  was  calculated  to 
give  1  atm  of  pressure  at  the  operating  temperature.  These  runs  were 
filled  with  380  mm  HC1  at  room  temperature.  The  HC1  pressure  was  ^2  atm 
for  each  run  at  maximum  temperature.  The  increase  in  HC1  pressure  was 
suggested  by  results  from  our  previous  runs,  in  which  the  amount  of 
material  transported  had  been  relatively  low.  An  inspection  of  the 
reaction  equations,  _  _ _ -  -  . 

ZnS  +  In2S3  +  8HC1  =  ZnCl2  +  2InCl3  +  4H2  +  2S2  —  Znln^  +  8HC1, 

indicated  that  the  reaction  is  pressure  dependent,  and  that  it  may  be 
possible  to  increase  the  amount  of  transported  material  by  increasing 
the  pressure.  This  was  later  verified.  The  material  was  submitted  for 
X-ray  analysis.  The  compound  CdIn2S^  was  identified  by  X-ray  powder 
diffraction  analysis.  The  pattern  is  shown  in  Figure  18. 

c.  Crystal  Growth  of  ZnIn2S^ 

Experiments  to  synthesize  ZnIn2S^  were  conducted  once  again 
using  iodine  as  the  transport  medium.  The  amount  of  reaction  and 
transport  increased  significantly  in  runs  lasting  a  few  days  as 
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Figure  18.  X-ray  diffraction  pattern  of  Cdln^S^. 
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compared  to  previous  runs  (using  HC1)  lasting  several  weeks.  The 
transported  material  contained  several  compounds  of  Zn^In^S^  of  varying 
stoichiometry;  all  were  rich  in  In^S^y  and  some  free  ZnS  was  identified 
as  well  by  X-ray  powder  diffraction  techniques.  In  subsequent  runs, 
the  iodine  transport  technique  combined  with  addition  of  excess  ln2Sg 
(n/2%)  yielded  high-quality  platelets  of  Znln^S^.  The  yellow  mica-like 
pla  telets  were  identified  to  be  Znln^S^  by  X-ray  diffraction.  Some  small 
(mm  edge)  crystals  were  also  obtained. 

2.  Synthesis  from  Solution 
a.  CdIn2Tei!i 

Primary  emphasis  was  placed  on  the  crystal  growth  of  single¬ 
crystal  Cdln^Te^,  a  compound  that  we  earlier  had  observed  as  having  a 
high  value  for  the  low-frequency  dielectric  constant  from  solution.  Using 
the  solid-liquid  relationships  shown  (Figure  19)  in  the  pseudo-binary 
phase  diagram  CdTe-In^Te^  as  published^  in  the  literature,  we  attempted 
to  grow  single-crystal  material.  Melt  compositions  were  selected  lying 
in  the  composition  range  covered  by  the  line  a-b.  In  this  region,  solid 
Cdln^Te^  (3  =  solid  Cdln^Te^)  is  in  equilibrium  with  liquid  (richer 
in  In2Te^)  in  the  temperature  range  785°C  to  702°C.  By  selecting  a 
starting  composition  close  to  "a,"  crystals  of  pure  CdIn2Te^  can  be 
grown  by  lowering  the  temperature.  Below  702°C,  however,  additional 
phases  begin  to  solidify  and  cause  a  multiphase  polycrystalling  region 
to  begin  to  grow,  thus  preventing  further  single-crystal  growth. 

Several  runs  were  made  utilizing  the  conditions  discussed  above  to 
grow  single-crystal  CdIn2Te^.  Fused  Mb-mm-diameter  quartz  crucibles 
with  conical  bottoms  were  used  to  nucleate  single-crystal  growth.  The 
sealed,  evacuated  ampoules  were  processed  in  a  vertical  two-zone 
Bridgman-like  furnace.  The  ampoule  was  rotated,  and  vibrational  stirring 
was  applied  during  the  growth  cycle.  Several  ingots  were  produced  in 
which  the  lower  (conical-shaped)  sections  were  single  phase  and,  in  one 
case,  single-crystal  CdIn2Te^;  the  upper  parts,  as  anticipated,  were 
multiphase  and  polycrystalline.  The  compound  CdIn2Te^  was  identified 
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TEMPERATURE 


by  X-ray  diffraction  analysis  (Figure  20).  Laue  photographs  showed 
several  slices  to  be  single  crystals. 

X-ray  microprobe  analysis  of  the  single-phase  material  showed 
insignificant  variations  in  composition  from  nose  to  tail  with  an 
average  composition: 

58%  Te 
29%  In 
13%  Cd. 

There  are  indications  from  the  analysis  that  the  composition  is 

slightly  rich  in  In^Te^  and  that  there  is  some  solid  solubility  between 

CdIn_Te.  and  In„Te0.  However,  our  results  indicate  that  the  solid 
2  4  2  3 

solubility  line  between  the  regions  is  more  vertical  than  that  reported 
in  Ref.  6  (dashed  line  in  Figure  19)  and  probably  closer  to  the  vertical 
(dotted)  line  in  Figure  19. 

The  single-crystal  section  (cone-shaped)  and  the  final  multiphase 
section  are  shown  in  Figure  21.  The  boundary  between  the  sections  is 
noted.  The  tip  of  the  cone  appears  as-grown  and  indicates  a  poor  start 
with  nucleation  getting  underway  on  an  apparent  gas  void,  resulting  in 
a  concave  interface  observed  during  crystallization.  That  the  crystal 
grew  single  under  these  conditions  (see  Figure  22,  an  SEM  photograph 
of  the  as-grown  tip)  is  confirmed  by  the  growth  patterns  shown  in 
Figure  22  as  well  as  by  Laue  photography  (as  previously  indicated). 

D.  SYNTHESIS  OF  A^B1^1  COMPOUNDS 

A  polycrystalline  ingot  of  Cu2GeS^  was  obtained  by  quendhing  a  melt 
of  that  composition.  X-ray  diffraction  indicated  the  presence  of  at 

least  one  additional  phase.  Table  3  shows  a  compilation  of  lattice 

1  IV  VI 

parameters  of  A_B  C_  compounds  with  literature  references  (compiled 

2  ^  3 

by  A.  Borshchevsky  of  Stanford  University)  and  indicates  the  existence 
of  solid  phase  transitions  and  a  variety  of  opinions  on  crystal  class 
of  the  low-temperature  form.  Both  the  ingot  and  the  starting  powder 
were  annealed.  After  annealing,  no  changes  were  observed  in  the  powder. 
Published  phase  diagrams  indicate  that  this  compound  is  a  peritectic  and 
can  be  grown  from  a  GeS2~rich  melt.  An  X-ray  diffraction  pattern  for 
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Figure  20.  X-ray  diffraction  pattern  of  CdIn0Te^, 


Figure  21.  Ingot  containing 
(Region  to  right 


graph  of  tip  of  Cdln  Te 


Table  3.  Lattice  Parameters  of  I^IV  VI^  Compounds 


Compound 

Crystal 

System 

a, 

X 

X’ 

B. 

deg 

Reference 

Cu2SiS3  (H.T.) 

Hex 

3.684 

6.044 

8 

Cu2SIS3  (L.T.) 

Tetr 

5.290 

5.078 

8 

Ortho 

11.21 

12.04 

6.03 

9 

Mono 

11.51 

5.34 

8.16 

98.95 

9 

Ortho 

10.981 

6.416 

6.046 

Cu2SiSe3 

Mono 

12.10 

5.62 

8.61 

99 

9 

Cu2S[Te3 

Mono 

12.86 

6.07 

9.05 

99 

9 

Cubic 

5.93 

8 

Cu2CeS3  (H.T.) 

Cubic 

5.  317 

8 

Cu2CoS3  (L.T.) 

Tetr 

5.219 

8 

Tetr 

10.41 

9 

Ortho 

11. 321 

3.766 

5.21 

10 

Cubic 

5.30 

11 

Mono 

6.433 

11.300 

7.533 

125.17 

12 

Mono 

7.464 

22.38 

10.64 

91.17 

13 

Cu2CeSe3 

Tetr 

5.595 

5.482 

8 

Tetr 

5.590 

10.97 

9 

Ortho 

5.591 

5.562 

5.488 

14 

Cubic 

5.55 

11 

Tetr 

5.591 

5.485 

15 

Ortho 

11.860 

3.960 

5.485 

10 

Cu2CeTe3 

Tetr 

5.956 

5.926 

8 

Tetr 

5.916 

11.85 

9 

Cubic 

5.95 

11 

Cu2SnS3 

Cubic 

5.445 

8 

Tetr 

5.426 

10.88 

9 

Cubic 

5.43 

11 

Cu2SnSe3 

Cubic 

5.696 

8 

Tetr 

5.689 

11.37 

9 

Cubic 

5.68 

11 

Cubic 

5.6877 

15 

Cu2SnTe  3 

Cubic 

6.047 

8 

Tetr 

6.048 

12.11 

9 

Cubic 

6.04 

11 

7209 
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Ci^GeS^  (Figure  23)  shows  the  X-ray  diffraction  pattern  of  the  high- 
temperature  distorted-cubic  phase  quenched  from  the  melt  (Figure  23(a)); 
the  annealed  powder  pattern  (450°C  for  500  hr)  is  shown  in  Figure  23(b). 

The  phase  transition  is  reported  to  be  at  670°C.  Although  there  is  a 
slight  difference  in  the  patterns,  we  have  not  observed  any  definite 
evidence  of  a  low-temperature  phase. 

The  structure  of  C^GeTe^  as  indicated  by  X-ray  diffraction  analysis 

O 

appears  to  be  cubic  with  a  =5.95  ±  0.02  A.  No  evidence  of  the  reported 
tetragonal  phase  is  indicated  in  the  X-ray  diffraction  patterns 
(Figure  24)  for  as-grown  and  annealed  (450°C  for  500  hr)  samples.  How¬ 
ever,  differential  thermal  analysis  (DTA)  indicates  the  existence  of  a 
solid  phase  transformation  (exact  temperature  undetermined)  and  a  con¬ 
gruent  melting  compound  at  490°C.  The  material  contains  a  small  amount 
of  a  second  phase  both  before  and  after  annealing.  Although  the  litera¬ 
ture  (see  Table  3)  is  not  precise  on  that  matter,  there  are  reports  of 
the  observation  of  a  tetragonal  phase.  Additional  DTA  work  is  required 
if  the  compound  appears  interesting  (i.e.,  has  a  high  dielectric 
constant) . 

E.  QUARTERNARY  COMPOUNDS:  Cu2CdCeTe4 

As-grown  and  annealed  samples  of  the  reported  quaternary  compounds 
(Cu^CdGeTe^)  both  show  by  X-ray  diffraction  analysis  (Figure  25)  free 
CdTe.  The  annealed  sample  (450°C  for  500  hr)  shows  (Figure  25(b))  reduced 
CdTe  peaks  and  a  shift  (change  in  lattice  parameters)  from  the  as-grown 
sample  (Figure  25(a)),  which  may  indicate  that  continued  annealing  will 
yield  the  quarternary  from  a  composition  that  looks  (as-synthesized. 

Figure  25(a))  like  the  ternary  Cu2GeTe2  plus  CdTe.  The  presence  of 
additional  phases  in  the  material  also  was  observed  both  before  and 
after  annealing;  perhaps  this  indicates  two  additional  phases  plus  CdTe. 
DTA  data  showed  at  least  three  phase  transformations:  786,  523,  and  513°C. 


(a)  As-grown  sample 


(b)  After  annealing 

Figure  23.  X-ray  diffraction  pattern  of  C^GeS^. 
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Figure  24.  X-ray  diffraction  pattern  of  cubic  (X^GeTe^. 
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(a)  As-grown  sample 


9001— S 


(b)  After  annealing 

Figure  25.  X-ray  diffraction  pattern  of  CX^CdGeTe^. 


48 


' 


M12680 


Figure  27.  Electrical  testing  apparatus  (internal  view). 
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SECTION  4 


MATERIALS  CHARACTERIZATION 

A.  MEASUREMENT  OF  DIELECTRIC  CONSTANT  AND  DIELECTRIC  RELAXATION 

The  ability  to  make  accurate  capacitance  and  dielectric  relaxation 
measurements  on  a  new  insulating  or  semi-insulating  compound  depends 
markedly  on  the  availability  of  a  good  metalization  system  for  electri¬ 
cal  contacts.  If  high  Schottky  barriers  are  formed  at  the  metal- 
contact/insulator  interface  and  the  insulator  does  not  have  a  suffi¬ 
ciently  high  resistance  (10^  Q  or  higher),  the  measured  capacitance  will 
be  rather  high  and  essentially  will  represent  the  Schottky-barrier  capa¬ 
citance.  This  phenomenon  was  observed  with  ZnSiAs^  having  silver  paint 
electrodes.  An  absurdly  high  value  (e  =  9500)  was  measured  for  the 
dielectric  constant.  That  this  was  genuinely  a  Schottky-barrier  effect 
was  ascertained  by  applying  a  voltage  bias  to  the  capacitor  and  observ¬ 
ing  the  change  in  capacitance.  Although,  for  any  given  material,  a  metal 
or  an  alloy  system  may  be  found  to  yield  non-barrier-forming  contacts, 
the  search  for  such  contacts  for  the  new  materials  made  during  this 
program  would  represent  a  major  undertaking.  Alternatively,  if  a  suit¬ 
able  metal  that  is  also  a  constituent  of  the  ternary  chalcogenide  (e.g., 
In  for  Cdln^Se^,  CuInSe^,  and  AglnSe^)  is  used  as  the  contact  and  fur¬ 
ther  diffused  into  the  shallow  surface  layers  of  the  chalcogenide,  a 
graded,  nearly  Ohmic  junction  may  be  expected  at  the  metal-insulator 
interface.  Although  this  approach  has  worked  well  with  single  crystals 
(e.g.,  Ag  on  AgGaS^),  it  has  one  main  drawback  with  polycrystalline 
insulators.  Even  after  a  moderate  diffusion  anneal  of  In  in  polycrystal 
line  CuInSe^,  AglnSe^,  and  Cdln^Se^,  we  observed  that  the  samples  were 
essentially  shorted  or  exhibited  resistances  on  the  order  of  a  few 
hundred  ohms  and  therefore  could  not  be  measured.  We  believe  that  the 
shorts  are  caused  by  diffusion  of  the  metal  along  grain  boundaries  and 
the  subsequent  formation  of  conducting  filaments  across  the  bulk  of  the 
insulator. 
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As  a  solution  to  these  problems,  we  decided  that  capacitance 

measurements  should  be  made  with  two  blocking  contacts.  This  also 

prevents  carrier  injection.  The  blocking  electrodes  consist  of  two  thin 

(0.004  in.)  sheets  of  Mylar  with  evaporated  silver  contacts.  The  Mylar 

capacitors  have  the  same  area  as  the  test  insulator.  The  insulator  is 

sandwiched  between  the  two  Mylar  capacitors,  and  the  capacitance  C  and 

m 

the  equivalent  parallel  resistance  R  of  the  sandwich  structure  at  a 

m 

frequency  10  are  measured  using  a  Boonton  750C  capacitance  bridge.  The 
insulator  is  then  removed  from  the  sandwich,  and  the  total  capacitance 
Cg  and  shunt  resistance  are  measured  at  the  same  frequency,  0).  A 
network  analysis  of  the  equivalent  circuits  for  the  two  cases  yields  the 
capacitance  and  resistance  R^  of  the  test  insulator  (as  shown  in 
Appendix  A. 

The  frequency  range  was  from  5  kHz  to  100  kHz.  These  measurements 
were  performed  at  room  temperature  using  peak  test  signal  levels  of  1  to 

3  V. 


The  static  dielectric  constant  e  was  obtained  from  the  calculated 

s 

capacitance  (Appendix  A)  at  any  given  frequency  by  using  the  formula: 


E 

s 


V 


8.85  A 


9 


where  d  is  the  thickness  of  the  capacitor  (in  m) ,  A  is  the  area  of 

2 

coincidence  of  the  electrodes  (in  m  ),  and  is  the  measured  capacitance 
(in  PF). 

A  general-purpose  portable  electrical  test  fixture  was  fabricated 
for  dielectric  constant,  dissipation,  resistivity  (using  Van  der  Pauw's 
four-point  probe  method),  dielectric  breakdown,  and  EO  coefficient. 
Measurements  are  possible  with  this  apparatus  at  any  temperature  from 
77°C  to  600°C.  In  addition,  the  fixture  can  be  evacuated  to  10  ^  Torr. 
Figures  26  and  27  show  some  of  the  salient  features  of  this  facility. 

The  measurements  made  on  materials  synthesized  during  this  program  are 
summarized  in  Table  4.  Data  on  dissipation  factor  (tan  6)  is  given 
in  Table  5. 
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Table  4.  Dielectric  Constants  at  Indicated  Frequencies  for  Some 
Binary  and  Ternary  Chalcogenides 


Compound 

Dielectric  Constant 

Frequency,  kHz 

GeS2 

10.7 

10 

GeSe  2 

5.32 

6 

5.31 

10 

5.30 

50 

5.30 

100 

5.26 

400 

AgGaS2 

13.7 

10 

AgGaTe2 

200 

5 

270 

1000 

ZnSiAs2 

34.78 

6 

34.65 

10 

33.64 

100 

32.91 

400 

ZnGa„S. 

2  4 

62.1 

6 

40.9 

10 

35.8 

20 

21.3 

50 

18.1 

100 

13.5 

400 

CuInSe2 

36.81 

6 

36.15 

10 

35.74 

100 

35.71 

400 

Cu2GeS^ 

49 

10 

Cu0CdGeTe , 

2  4 

60.3 

10 

Cdln  Te . 

2  4 

264.5 

10 

7209 


Table  5.  Dissipation  Factor  (tan  5) 
for  Several  Chalcogenide  Materials 


Material 

tan  6 

ZnGa0S. 

2  4 

0.499 

Cu2GeS3 

0.00028 

Cu2CdCeTe^ 

0.00364 

Cdln^Te . 

2  4 

0.00515 

mm 


Numerous  samples  were  synthesized  that  could  not  be  measured 
because  their  resistivities  were  too  low  for  our  set-up.  These 
included: 


AgInSe2 

(Cu,Ag)GaS2 

AgGa (S , Se) ^ 

AgInTe2 

TlSe 

CdIn_S, 

2  4 

Ag(Ga,In)S2 

Ag^AsS^ 

•  AnIn„Se. 

2  4 

Unfortunately,  the  scope  of  the  program  did  not  allow  detailed 
study  of  these  materials  to  allow  us  to  control  resistivity  (e.g., 
annealing,  doping  studies).  This  is  discussed  in  more  detail  in 
Section  6.  In  addition,  absorption  caused  by  free  carriers  precluded 
far-IR  spectroscopy  of  many  of  these  samples. 

B.  DISCUSSION  OF  RESULTS 

Reproducibility  of  results  was  checked  by  measuring  two  separate 
samples  from  the  same  as-grown  ingot  at  two  different  low  frequencies. 


54 


The  reproducilbility  was  fairly  good.  The  small  sample-to-sample 
variations  noted  were  probably  caused  by  minor  local  variations  in 
structure,  impurity,  or  defect  concentrations. 

1.  GeS2 

Figure  28  shows  the  dielectric  constant  as  a  function  of  frequency 
in  GeS2*  The  change  in  the  dielectric  constant  is  fairly  small  (10.73  at 
5  kHz  as  compared  to  10.52  at  100  kHz).  Figure  29  shows  the  dissipa¬ 
tion  factor  for  GeS^  as  a  function  of  test  frequency.  The  curve  has 
some  interesting  features.  Two  strong  primary  peaks  are  seen  at  5.5 
and  65.  kHz;  these  are  followed  by  a  definite,  albeit  weaker,  secondary 
peak  at  45  kHz.  These  dielectric  relaxation  peaks  suggest  separate, 
distinct  relaxation  processes  and  deserve  further  investigation. 

The  data  shown  in  Figures  26  and  27  represents  measurements  on  only 
one  sample.  Since  the  as-grown  GeS^  crystal  boule  had  a  small  diameter 
and  the  crystal  tends  to  cleave  easily,  we  were  not  able  to  obtain  a 
large  enough  sample  on  which  large-area  electrodes  could  be  formed.  The 
small-electrode-area  sample  discussed  above  had  only  a  small  capaci¬ 
tance  (a  few  pF).  Because  the  capacitance  of  the  test  capacitor  was 
so  low,  accurate  measurements  of  the  dielectric  constant  and  dissipa¬ 
tion  factor  were  very  difficult.  Larger  crystals  which  would  yield 
larger  capacitance  (<100  pF)  capacitors  should  be  used  to  check  the 
reproducibility  and  accuracy  of  the  dielectric  constant  and  the 
relaxation  measurements. 

2,  GeSe2 

Figure  30  shows  the  variation  of  dissipation  factor  versus  fre¬ 
quency  of  GeSe^;  the  same  plot  for  ZnGa2S^  is  shown  in  Figure  31.  In 
the  case  of  GeSe2,  two  dielectric  relaxation  peaks  corresponding  to 
8.5  kHz  and  70  kHz  are  observed.  The  plot  for  ZnGa^S^  shows  only  a 
small  peak  at  27  kHz.  The  dielectric  constant  of  ZnGa^S^  is  much 
higher  at  low  frequencies  (6  kHz)  than  similar  constants  for  other 
materials  that  we  have  previously  investigated  in  this  project. 
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Figure  29. 


Dissipation  factor  versus  frequency.  Sample  is  GeS2* 
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Figure  30.  Tan  6  versus  frequency  for  GeSe^. 
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The  dissipation  factors,  on  the  order  of  1C>\  are  fairly  high  as 

compared  with  the  chalcogenides  (listed  above),  which  are  on  the 
_3 

order  of  10  .  At  present  we  do  not  understand  these  relaxation  phenom¬ 

ena  because  of  a  lack  of  knowledge  of  the  details  of  the  local  atomic 
arrangements  and  defect  structures  in  this  material. 

3.  AgGaTe2 

The  0. 5-in. -diameter  boule,  grown  by  the  Bridgman-Stockbarger  tech¬ 
nique,  was  found  to  be  polycrystalline  along  its  entire  length.  A  few 
large  pieces  were  mined  out  for  dielectric  measurements.  The  average 
dielectric  constant,  measured  at  1  MHz  for  this  material,  was  270 
(Table  4).  However,  there  were  large  variations  in  the  measurements 
from  sample  to  sample  (nearly  35%),  indicating  variations  in  stoichiom¬ 
etry.  IR  transmission  measurements  (Figure  32)  showed  transmission  of 
^40%  (uncorrected  for  Fresnel  losses)  at  2.5  ym,  decreasing  in  a 
manner  typical  of  free-carrier  absorption,  becoming  opaque  beyond  4.5  ym. 
From  its  composition  and  relationship  to  other  chalcopyrites ,  trans¬ 
mission  further  out  in  the  IR  would  be  expected.  DC  resistivity  for 
AgGaTe2  was  measured  to  be  ^6  x  10^  fl-cm. 

4.  CdIn2Te^ 

The  dielectric  constant  measurements  and  relaxation  spectra  for 
this  material  are  shown  in  Tables  4  and  5.  We  attempted  to  detect 
f erroelectricity  in  this  material  by  employing  thermally  induced 
depolarization  currents  in  the  range  of  77  to  300°K.  This  work  was 
done  under  the  direction  of  Professor  Ruben  Braunstein,  Department 
of  Physics,  UCLA.  No  peaks  in  the  depolarization  current  versus 
temperature  were  observed,  indicating  the  absence  of  the  effect  at 
least  in  the  above  temperature  range.  Similar  measurements  at  high 
temperatures  (300  to  750°K)  also  showed  no  peaks. 
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The  dielectric  relaxation  spectrum,  a  plot  of  the  dissipation 
factor  (tan  <5)  versus  frequency  for  ZnSiAs^,  is  shown  in  Figure  33. 

There  are  two  distinct  peaks  corresponding  t6  frequencies  of  60  kHz  and 

k 

220  kHz,  which  correspond  to  relaxation  times  of  16  psec  and  0.7  usee, 
respectively.  Although  sufficient  details  of  the  defect  structure  of 
this  sample  are  not  presently  available  to  draw  any  definite  conclu¬ 
sions,  we  can  speculate  that  the  shorter  time  constant,  0.7  ysec, 
corresponds  to  the  jump  frequency  of  the  lightest  of  host  ions,  silicon, 
into  and  out  of  Si  vacancies.  The  longer  time  constant  of  16  ysec  may 
be  due  to  similar  jump  frequencies  of  the  other  host  cation,  zinc.  The 
host  anion,  arsenic,  because  of  its  large  mass,  is  not  expected  to 
respond  to  the  applied  high  frequency  of  electric  fields  corresponding 
to  the  observed  peaks. 

« 

C.  FAR-IR  SPECTROSCOPY 

Far-IR  spectroscopy  in  the  range  of  5  to  100  cm  ^  was  performed  on 

all  synthesized  materials  that  were  transparent  in  this  section  at 

these  frequencies.  Transmission  spectra  were  obtained  and  analyzed 

under  the  direction  of  Professor  Paul  L.  Richards,  Department  of 

Physics,  University  of  California  at  Berkeley.  The  index  of  refraction 

is  estimated  from  the  measured  mechanical  thickness,  t,  and  the  observed 

wavelengths  of  the  maximum,  A  ,  and  minimum,  A  .  ,  in  transmission: 

max  min 

\  «=  2nAt 

max  m 

\  =  4nAt 

min  2m+l  * 

where  m  is  any  integer  21.  These  equations  are  applied  in  the  spectral 
regions  where  the  complex  part  of  the  index  is  negligible.  At  these 
frequencies,  the  dielectric  constant,  c,  is  equal  to  the  index  of 
refraction  squared: 

2 

e  «  n 
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Figure  33.  Dielectric  relaxation  of  ZnSiAs2» 


1.  GeS2 

The  far-IR  spectrum  of  GeS^  is  shown  in  Figure  34.  The  top  curve 
(1)  in  Figure  34  shows  the  transmittance  through  a  GeS2  sample  0.356  mm 
thick  at  a  temperature  of  1.3°K.  The  spectrum  shows  interference 
fringes  and  sharp  absorptions  at  48  to  63  cm  \  Curve  2  shows  the 
estimated  index  of  refraction.  The  calculated  dielectric  constant  from 

this  data  is  in  good  agreement  with  measured  values.  The  dispersion  due 

to  the  48  to  63  cm  ^  modes  and  also  to  higher  modes  shows  clearly  in  the 

diagram. 
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2.  AgInTe2 

The  far-IR  transmission  spectrum  of  AglnTe^  (Figure  35)  was 
determined  and  analyzed  from  a  sample  with  a  measured  thickness  of 
0.113  cm.  The  index  of  refraction  was  determined  to  be  3.40,  and  no 
noticeable  dispersion  over  the  measured  frequency  range  was  observed. 

3.  Cdln^e^ 

Far-IR  transmission  measurements  on  Cdln  Te,,  without  AR  coatings, 

-1  L  * 

in  the  range  of  500  to  40  cm  wave  numbers  (corresponding  to  20  to 
2500  ym)  show  a  constant  relative  transmission  of  20%  between  20  to 
32  ym,  followed  by  a  medium  absorption  peak  at  340  cm  ^  and  two  strong 
absorptions  at  212  cm  ^  and  at  125  cm  ^  (Figure  36).  The  nature  of 
these  absorptions  is  not  yet  understood  and  will  require  further  trans¬ 
mission  studies  at  low  temperatures  with  increased  instrument  resolution. 

4.  AgInSe2 

The  spectrum  for  a  29.4-mil-thick  sample  of  AglnSe^  is  shown  in 
Figure  37.  An  average  index  of  3.18  in  the  indicated  range  was  deter¬ 
mined  by  this  measurement. 

5.  TISe 

TISe  of  18.0-mil  thickness  (Figure  38)  was  determined  to  have  an 
average  index  of  4.63  in  the  range  15  to  23  cm 

6.  Mixed  Crystal  AgGa(S,Se)9 

A  sample  of  the  mixed  crystal  AgGa)S,Se)2>  where  [S]  =  [Se],  having 

a  thickness  of  58.9  mil  was  also  evaluated  in  the  spectral  region  of  0  to 

40  cm  1  (Figure  39).  The  index  of  refraction  between  7  and  26  cm  ^  was 

determined  to  be  2.93  ±  0.05.  Measured  EO  coefficient  values  of  the 

-2 

end  members,  AgGaS„  and  AgGaSe„,  have  been  reported  at  3.2  x  10  M/V 

-12  1  1 
and  3.8  x  10  M/V,  respectively. 
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TRANSMITTANCE 


9324-1 


Figure  36.  Far-IR  transmission  of  Cdlr^Te^. 
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.  Far-IR  transmission  of  AgInSe2» 


Figure  37 
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Figure  38.  Far-IR  transmission  spectrum  of  TISe. 


7 


Other  Samples 


Samples  of  CuInSe^  and  Ag^AsS^  with  thicknesses  of  0.016  and 

0.014  in.,  respectively,  showed  no  far-IR  transmittance  at  a  sensitivity 

-3  -1 

level  of  T  <  10  over  the  frequency  range  10  <  V  <  40  cm  . 

All  that  we  can  conclude  is  that  there  is  either  very  large  absorp¬ 
tion  in  these  samples  or  that  the  index  of  refraction  is  much  greater 
than  10,  giving  very  high-0  Fabry-Perot  resonances  that  are  smeared  out 
by  irregularities  in  the  sample  thickness. 

Additional  samples  submitted  for  measurement  were  not  transparent 
in  the  far  IR,  Several  low-resistivity  samples  with  apparent  free- 
carrier  absorptions  could  not  be  analyzed  by  our  methods. 
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9524—6 


Figure  39.  Far-IR  transmission  of  AgGa(S,Se)2«  Reference  spectrum  is 
scan  without  sample.  Peaks  above  and  below  reference  level 
indicate  channel  spectra  due  to  resonances  in  sample. 

Ratio  represents  transmittance  (i.e.,  ratio  of  sample  spec¬ 
trum  (point-by-point)  to  reference  spectrum). 
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SECTION  5 


THEORETICAL  MODEL:  ELECTROOPTIC  EFFECT  IN  CRYSTALS 


During  this  program,  a  parallel  effort  was  pursued  by 
Professor  Amnon  Yariv  and  Dr.  C.  Shih,  Applied  Physics  Department  at 
Caltech,  to  develop  a  theroetical  model  for  the  EO  effect  in  coordina¬ 
tion  with  the  structural  analog  concepts  pursued  in  this  program.  The 
prime  motivation  of  the  theoretical  effort  was  to  obtain  an  expression 
that  could  be  used  to  calculate  and  predict  the  EO  coefficient  of  new 

crystals.  Yariv  and  Shih,  by  applying  the  localized  bond  charge  model 

17  18 

o.f  Phillips  and  Van  Veehten,  *  which  attributes  the  dielectric 
response  of  covalent  crystals  to  the  localized  bond  charge,  have 
developed  a  method  for  calculating  the  EO  tensor  coefficients.  Compari¬ 
sons  with  experimental  values  for  binary  compounds  having  zincblende 

t 

and  wurtzite  structures  were  excellent.  Their  method  and  results  are  ; 

presented  in  detail  in  Appendix  B  (Ref.  19).  i 

i 

Yariv  and  Shih  recently  extended  the  calculation  to  ternary  com¬ 
pounds  (e.g.,  LiNbO^ ,  Cdln^Te),  where  they  also  find  very  good  agree¬ 
ment.  These  results  are  presented  in  Table  6. 

.We  applied  the  theoretical  model  for  calculating  the  EO  coeffi- 

i  . 

cient  to  the  crystal  Cdln^Te^.  We  used  the  following  data: 

a  'Covalent  radii  (Cd.In.Te)  —  1.405  A  . 

J  J  1 


Structure  —  4  - 


Using  published  data  of  bond  susceptibilities,  we  calculated  the 
optical  dielectric  constant  =  7.25.  Using  the  above  data  in  the 
equations 


(r  ) 

ijx  ionic 


e  (£*  .  -  e  . ) 

o  dck  c°k 

*  '  • 

V  Ne  £.£, 
c  i  j 


I  £ 


Table  6.  Results  for  LiNbO^  and  LiTaO^ 


exptl 


LiNb03 

LiTa03 

2t 

4: 

4; 

43 

3 

-.8 

/■ 

2 

!.0 

Nb-0  (Short) 

Nb-0  (Long) 

Ta-0  (Short) 

Ta-0  (Long) 

1.889  A 

2.112 

1.891 

2.071 

0.821 

0.830 

0.847 

0.853 

0.292 

0.241 

0.282 

0.238 

r33 

r51 

r33 

r51 

+19.9 

+19.7 

+27.8 

+16.5 

+  6.0 

+  0.8 

+  3.70 

+  0.24 

+25.9 

+20.5 

+31.5 

+16.7 

+28 

+23 

+30 

+15 

-  1.48/  fj,  -  0.02  , 


we  obtained 


220  x  10-12  ra/V  ±  30%  . 


This  predicted  value  is  a  huge  number  (roughly  100X  that  of  GaAs  and 
10X  that  of  LiNbO^).  The  main  difference  between  this  crystal  and 
LiNbO^,  as  an  example,  is  due  to  the  factor: 
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eiej 


which  is  equal  to  8.55  in  Cdln^Te.  and  to  1.02  in  LiNbO_.  It  will  be 

2  k  3  - 

very  exciting  to  see  if  this  prediction  is  actually  true. 

We  also  predict  a  large  EO  coefficient  for  Anln^S^  since  it  has 
3ra  symmetry  so  that  the  form 


l  -;4jk 


u 


is  non-vanishing  for  r^  .  The  more  numerical  prediction  awaits  the 

result  of  measurement  of  e,  . 
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SECTION  6 


CONCLUSIONS  AND  RECOMMENDATIONS 


At  the  conclusion  of  this  study  program,  there  are  many  observations 
and  experiences  to  be.  noted.  We  have  had  a  very  successful  program  which 
has  illuminated  several  possible  materials  having  potentially  high  EO 
coefficients.  These  materials,  especially  Cdln^Te^  and  AgGaTe^,'  have 
been  shown  to  have  high  dielectric  constants,  as  we  had  predicted  from 
the  structural  analog  model.  In  addition,  a  new  theoretical  calcula¬ 
tion  (Appendix  A)  is  now  available  for  the  calculation  of  the  EO  coeffi¬ 
cients  based  on  the  localized  band  charge  model  of  Phillips  and  Van  Vechten. 
This  was  developed  in  a  parallel  effort  by  A.  Yariv  and  C.  Shih  of  Caltech. 

Many  problems  arose  during  this  program  concerning  materials 
synthesis,  crystal  growth,  and  measurements.  We  were  grossly  misled  by 
many  reports  inthe  Russian  literature,  especially  concerning  crystal 
synthesis  by  vapor  transport  techniques.  Fortunately,  our  association 
with  Dr.  Alexander  Borshchevsky,  formerly  at  the  Center  for  Materials 
Research,  Stanford  University,  led  to  the  clarification  of  the  Russian 
literature  both  as  to  content  and  quality. 

In  addiition  to  misleading  literature,  the  ternary  chalcogenides 
have  some  unique  synthesis  problems.  Particularly  difficult  are  combina¬ 
tions  where  one  of  the  binary  constituents  (e.g.,  ^n2S3  or  *s 

considerably  more  volatile  than  the  other  (e.g.,  ZnS  or  CdS).  This 
typically  creates  problems  in  achieving  the  desired  stoichiometry. 

An  extreme  example  is  Znln^S^,  where  we  encountered  approximately  a  dozen 
identified  compounds  of  the  type  Zn^In^S^.  The  desired  stoichiometry, 
Znln^S^,  was  finally  achieved  by  regrinding  the  other  compounds  and 
adding  excesses  of  and  S2. 

There  are  many  areas  of  study  (pointed  out  above)  where  considerably 
more  work  is  needed.  Surveying  a  wide  range  of  compounds,  especially 
"new"  compounds  where  there  is  little  information  on  synthesis  or  on 
properties,  only  serves  to  open  more  areas  where  research  is  needed. 

We  have  only  scratched  the  surface  to  families  of  new  materials  with 
a  wide  variety  of  applications. 
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Many  crystalline  samples  were  fabricated  but  could  not  be  measured 
because  of  their  low  resistivity.  These  samples,  in  addition,  exhibited 
new  transmission  in  the  far  IR  due  to  free  carrier  absorption.  Some 
annealing  was  attempted  to  increase  resistivity  by  modifying  crystal 
stoichiometry,  but  the  scope  of  this  program  would  not  allow  an  in 
depth  study  of  these  phenomena. 

The  high  dielectric  constant  and  structural  analog  approach  which 
led  to  the  selection  of  new  materials  with  potentially  higher  EO  coeffi¬ 
cients  also  carries  with  it  the  possibilities  of  higher  nonlinear  coeffi¬ 
cients  and  higher  acoustooptic  figures  of  merit  compared  to  known 
materials.  The  applications  of  these  materials  go  beyond  that  which 
was  our  main  objective  —  the  EO  tunable  filter.  New  materials  in  these 
families  can  be  selected  with  wide  variations  in  band  gap  and  phase 
matching  from  which  to  select  new  improved  materials  for  frequency  upcon- 
version,  parameteric  oscillators,  SAW  devices,  modulators,  and  IR 
detectors. 

We  recommend  that  the  study  of  these  ternary  compounds  be  con¬ 
tinued  with  the  applications  in  mind  so  that  standard  synthesis  and 
crystal  growth  techniques  can  be  developed  to  make  them  available  for 
future  devices. 

In  many  aspects  of  the  program,  we  have  made  measurements  and 
observations,  but,  because  of  the  limited  depth  of  study  required  in  any 
program  of  limited  resources,  we  were  unable  to  draw  significant 
conclusions.  One  prime  example  is  the  work  done  in  dielectric  relaxa¬ 
tion,  where  an  in  depth  study  of  defects  of  the  crystal  would  be 
required  to  understand  and  interpret  the  data. 

In  conclusion,  we  have  developed  a  new  approach  which  is  now 
backed  by  theory  for  the  selection  of  new  materials  with  potential  for 
higher  performance  in  a  variety  of  devices.  This  field  should  now  be 
pursued  in  greater  depth. 
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APPENDIX  A 

ANALYSIS  OF  CAPACITANCE  MEASUREMENTS  USING  BLOCKING  CONTACTS 


Problem; 

An  impedance  bridge  measures  the 
total  capacitance  CM  and  the  total 
resistance  (equivalent  shunt)  of 
the  two  capacitors  C5  and  C^  in  series 
at  a  frequency  f.  The  bridge  inde¬ 
pendently  also  measures  Cg  and  Rg. 

Find  the  value  of  Cj  and  R^. 


8463-4 


Ps 

Ct 


Solution ; 


Since  Cg  and  Rg  are  in  parallel,  the  net  (total)  impedance  Zg  of 
Cg  and  Rg  is  given  by 


_L  l_L  +  V) 
zs  I  Rs  ^  I' 


where  to  =  2iTf,  and  j  =  /^f.  Therefore,  it  follows  that 


2  2  2 

d  +  m  RgCg) 


RS  V 


2  2  2 

d  +  .  RsCg) 


Similarly  for  the  circuit  of  C^,  and  R^, 


+  j 


CTm 


T  2  2  2  '  J  222 

1  +  I  +  c/r^ 


Since  Z  and  Z  are  in  series,  the  total  measured  Zw  is  equal  to 
o  i  M 

Zc  +  Z_.  The  measurement  of  Z,,  yields  a  value  for  R,,  and  Cw  since 
o  i  n  MM 


'M 


-j— t  +M 

I  "m  j  i 
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j  trrl  l 


T 


+  j 


2  2  2 
^  Vh 


,  ,  2  2  2 
14w  Vm 


From  =  Zg  +  Z^,,  we  get,  after  equating  real  and  imaginary  parts. 


*M 


R„ 


2  2  2 
1  **»  Vm 


2  2  2 
i-h,  RSCS 


+ 


*T 


2  2  2 
1-Hj  RC~ 
T  T 


(A-l) 


and 


j 

I* 


\ 

■I 


^CM 


ESCS 


2  2  2 

1+“  Vm 


2  2  2 
i-Ko  RSCS 


+ 


4ct 


2  2  2 
14u  RXC^ 


Therefore , 


(A-2) 


*T 


*M 


2  2  2 
1-h/rjcJ 


2  2  2  2  2  2 
1-Ho  R>1Cm  1-H,  RsCs 


=  8,  a  constant. 


since  R^,  R^,  C^,  C^,  and  ware  measured  values.  Then  it  follows  that 


2  2  2 

Rx  =  e  +  m  R^c^e  . 


Similarly,  from  Eq.  A-2  we  get 


r?ct 


RMCM 


r2scs 


2  2  2 

1+rn  R^ 


2  2  2 

1+w  Vm 


2  2  2 

1+w  Vs 


=  A,  a  constant 


(A- 3) 
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and 


R^CX  =  A  +  Aw2R^C2  .  (A -4) 

Thus,  for  the  two  unknowns  and  R^,  we  have  two  equations: 

R^c2w2e  -  rt  +  6  =  0  (A- 5) 

and 

R^Ct  -  X  -  Xu)2R^C2  =  0  .  (A- 6) 


The  simplest  solution  will  be  derived  from  Eq.  A-5: 


Rt  -  6 

~2  2~ 

Rp  w  e 


or 


Since  only  positive  roots  are  possible,  we  can  write 


„  1 
T  Rtw 

2 

Substituting  for  in  Eq.  A-5,  after  squaring  both  sides  of  Eq.  A-6 , 
yields 
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4  2  2  „ , 2  222  ,  2  4  4  4 

i^C^,  =  A  +  2A  u  RtCt  +  A  w  R^Ct 


From  this  we  obtain 


(Ve, 

^  R^2e 


-  ,  ,  R?  (Rt-0)  aVrJ(Rt-0)2 

Az  +  2AV  1 


2  2 
R^oi  0 


„4  4.2 
RTa)  0 


or 


2  ^  2  2A2 

— aT  *  x  +  —  <V9>  + 

0)  0 


X2(Rt-0)2 

7“" 


2  22222  22  2 
0R^(Rt  -0)  -  A  u  0*  +  2Xu  0(Rt-0)  +  Aw  (Rt-0) 


or 


229  9  92  9  99  22  2  22 

A  w  R^  -  2A  Rt0  +  re  +  2A  0Rt  -  2A  0  +  A  0  =  R^0-  R^.0 


Therefore , 


3  2  2  2  2  2 

R^0  -  Rt0  »  Xu  R^ 


giving 


2  2  2 

RT0  -  Aw  +0  , 


.2  2 

rt  -  +  0, 
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Quantitative  Calculation  of  Electro-Optic  Coefficients  of  Diatomic  Crystals 

Chun-Ching  Shih  and  Amnon  Yariv(l) 

Department  of  Applied  Pfiysics.  California  Institute  of  Technology,  Pasadena,  California  91125 

(Received  19  September  1979) 

The  bond-charge  dielectric  theory  of  Phillips  and  Van  Vechten  Is  applied  to  the  calcu¬ 
lation  of  the  electro-optic  tensor  coefficients.  The  agreement  of  the  theoretical  predic¬ 
tions  with  experimental  values  in  tho  case  of  zinc  blende  and  wurtzlte  crystals  Is  very  . 
good. 


The  second-order  nonlinear  optical  response  of 
asymmetric  crystals  is  usually  represented  by 
the  relation  +  fi)=x  between  the  am¬ 

plitude  of  the  induced  polarization  at  u>  +  f7  and 
the  inducing  field  amplitudes  at  co  and  17.  The 
case  when  both  w  and  17  are  optical  frequencies, 
that  is,  frequencies  above  the  lattice  response 
but  below  optical  absorption,  has  been  considered 
by  Levine.1  He  used  the  localized  bond-charge 
model  [see  Fig.  1(a)]  of  Phillips  and  Van  Vechten2’3 
(PV)  which  attributes  the  dielectric  response  of 
covalent  crystals  to  the  localized  bond  charge  re¬ 
sulting  in  a  linear  susceptibility: 

X  =  (hSlp)2/Et2.  (1) 

Here  4lp  is  the  plasma  frequency  due  to  valence 
electrons  and  E f,  the  effective  energy  gap,  is  giv¬ 
en  by  E2  =  E2  +  C2,  where  E„  is  the  homopolar 
component  and  C  the  heteropolar  (ionic)  compo¬ 
nent  of  the  gap  energy.  Levine  starts  with  the 
linear  dielectric  response  P ^ =  XijEjw) ,  taking 


Td 

- 

- 

T-O 

l-j- 

1 

-< 

fa 

y-1-  © 

E  *  0  E  *  0 


(0)  (b)  (c) 


Cd) 


FIG.  1.  The  response  of  the  Ions  and  the  bond  charge 
to  the  applied  electric  field.  In  (a)  and  (b)  the  frequency 
of  the  field  Is  higher  than  the  lattice  rcsponso.  Only 
the  bond  charge  Is  displaced,  while  tho  tons  remain 
stationary.  In  (c)  and  (d)  the  frequency  of  the  field  Is 
lower  than  the  lattice  response.  It  Induces  a  displace¬ 
ment  of  tho  bond  charge  (A rn),  a  variation  of  the  bond  87 
length  (Adj),  and  a  bond  rotation  fid). 


Xu  to  be  an  instantaneous  function  of  the  second 
field  This  field  causes  a  change  Ara(=-Arfi) 

in  the  bond  charge  position,  as  shown  in  Fig.  1(b), 
which  oscillates  at  17.  The  explicit  dependence  of 
Et  and  C  on  ra  given  by  PV  is  then  used  to  obtain 
XtJ>,  where 

Xij  =  Xu0)+*X  (/(f) 

=x,/0)  +  2x„A(n)cosnt.  (2) 

If  the  frequency  (7  is  below  the  lattice  response 
region  (we  will  refer  to  it  in  this  case  as  “low”), 
then  in  addition  to  the  purely  electronic  nonlinear 
response  described  above  we  now  have  a  contri¬ 
bution  to  Xu*  due  to  the  fact  that  now  the  crystal 
ions  are  capable  of  following  the  field  £„<n>  cosfif. 
This  is  illustrated  in  Figs.  1(c)  and  1(d).  In  addi¬ 
tion  to  the  displacement  A ra  of  the  covalent  bond 
charge  we  now  have  an  elongation  A d0  of  the  atom¬ 
ic  separation  as  well  as  a  rotation  Ad  of  the  bond 
direction,  both  caused  by  the  ionic  displacement 
Ax*.  Ax*  is  obtained  from  "low”-frequency  dielec¬ 
tric  constant  measurements  and  is  used  to  deter¬ 
mine  d0  and  A 9.  We  use,  in  the  spirit  of  Levine,1 
the  change  A d0  to  calculate  the  corresponding 
change  Aycostlf  in  bond  susceptibility.  This  will 
give  rise  to  a  polarization  *  0)=XuajE' 

A  second  contribution  to  Xu»  is  due  to  the  rocking 
at  fl  of  the  bond  angle  (9  =  0o  +A  9  cosftf )  which 
yields  a  dipole  component  along  the  direction  i 
at  (w  +(7)  even  when  Ad0  =  0. 

In  what  follows  we  will  obtain  expressions  for 
the  ionic  contributions  to  xllk  which  are  due  to 
A d0  and  Ad.  When  the  result  is  added  algebraical¬ 
ly  to  Xu*el'c>  35  measured  by  second-harmonic- 
generation  experiments  or  calculated  by  Levine,1 
the  result  is  the  total  nonlinear  tensor  xUi 
=  Xuk°nic  +  Xu*clec  •  The  constants  Xu*  thus  de¬ 
termined  are  those  which  characterize  the  linear 
electro-optic  (Pockels)  effect.  The  relationship 
between  the  conventionally  defined  electro-optic 
tensor  rUk  and  x,j*  is 

Xijk =  ~  (c  (3) 

The  linear-  susceptibility  of  a  diatomic  crystal 

281 


©  1980  The  American  Physical  Socicfv 


Volume  44,  Number  4 


PHYSICAL  REVIEW  LETTERS 


28  January  1980 


is  given  as  x  =  {hQ,p)2/E2.2  The  ionicity  and  cova¬ 
lency  of  the  bond  are  defined  as  fi  =  C2/E2,  f  c 
^E2/E2?  The  expressions  used  in  the  evalua¬ 
tion  of  E k  and  C  are1'3 


Ekccro‘‘,  s  =  2.48,  (4) 

Cccexp(-kf0)(Za/ra-Z&/r&)e2,  (5) 

ra^r^r0  =  d0/2, 


where  d0  =  ra+r&  is  the  bond  length,  ya  8  are  the 
atomic  radii,  and  exp(-&/0)  is  the  Thomas-Fer- 
mi  screening  factor.  In  order  to  consider  crys¬ 
tals  with  highly  unequal  atomic  radii,  a  general¬ 
ized  form  for  E A  was  proposed1  as 


Ek-*r 


„[(ra-rr)23  +(r s  -rr)2 
’  2(r0-re)23 


(6) 


where  r e  is  the  average  core  radius. 

The  lineai-  macroscopic  susceptibility  tensor 
Xu  is  related  to  the  bond  polarizability  0„  by 

XiJ=V'tLnCtnianA,  (?) 


where  V  is  the  volume  of  unit  cell,  ani  is  the  di¬ 
rection  cosine  of  the  nth  bond,  and  the  summation 
is  over  all  the  bonds  in  a  unit  cell.  Although  PV 
describe  the  macroscopic  susceptibility  x  in 
terms  of  the  average  energy  gap,  we  assume  that 
and  C  are  also  related  to  the  bond  polarizabili¬ 
ty  directly,  i.e. , 

0(,cc  (hSlp)2/E  2 .  (8) 

When  the  bond  length  varies,  it  is  reasonable  to 
assume  that  the  ratio  of  ra  and  r6  remains  con¬ 
stant.  With  this  assumption,  the  two  independent 
parameters,  ro  6,  can  be  transformed  into  two 
quantities  which  relate  directly  to  the  macroscop¬ 
ic  properties  of  crystals: 

=(ra/d0)Ad0  +  6 

Ar6  =  {rB/d0^d0-i,  (9) 

where  6  is  the  displacement  of  the  bond  charge  in 
the  case  of  no  bond  elongation,  Ado  =  0,  From 
Eqs.  (5),  (6),  (8),  and  (9),  the  change  of  bond  po¬ 
larizability  is  obtained  as  (here  we  drop  the  bond 
index  n) 


0 


Zc,  +  Z  8 
Za  —  Z  g 


+  .S  (2s 


-1) 


frPd 

(>*0  -r 


2 

SL. 


c> 


4 


(10) 


where  p  =  (ra -r  6)/(ra +r  A  In  the  first  term,  hf0/2  is  obtained  because  the  screening  wave  number 
kt  is  proportional  to  d„~l  ,  and  the  number  -1  is  because  is  proportional  to  d0"j/2.  The  second 
term  on  the  right-hand  side  of  (10)  is  identical  to  that  obtained  by  Levine  in  his  calculation  of  the  non¬ 
linear  optical  susceptibility.5  The  first  term,  which  is  proportional  to  Atf0,  is  thus  the  ionic  contribu¬ 
tion  of  a  single  bond  due  to  bond  stretching. 

The  rotational  contribution  can  be  obtained  by  considering  the  changes  in  bond  direction  cosines. 
These  are  related  to  the  ionic  displacement  Ax,  by  Aan(  =  (6u- 


From  (7)  we  have 

AXu=V'lSn(an(QuA0n+Aan(a„A+an(AQni0j-  (H) 

The  complete  ionic  contribution  to  the  nonlinear  susceptibility  is  thus 

AXu'°n  =  {S-((5-/v,roX/a1,(Q„ja„*  +  i(aB(6j*  +  aBj64J}A^,  (12) 

where 

/“/ill  +  i(*,roS  +«/,  -  2.5  =[(U/0)  -  1.48]/,  -  0.02.  (13) 

Axt  is  related  to  the  dielectric  constant  of  the  crystal  as 

Wc/Ax*  =  £„(£</-£„')£„  (14) 


where  A'  is  the  number  of  pairs  of  atoms  per  unit  cell,  ec *  is  the  Callen  effective  ionic  charge*  edc'  is 
the  relative  dielectric  constant,  cj  is  the  relative  optical  permittivity,  and  E„  is  the  low-frequency 
electric-field  component  along  the  k  direction. 

Using  (3)  we  obtain  the  final  working  expression  for  the  Ionic  component  of  the  electro-optic  tensor: 
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TABLE  I.  Comparison  of  calculated  and  measured  electro-optic  coefficients  - - 

(clamped)  of  zinc  blende  (>  14)  and  wurtzite  (r33)  crystals  (In  units  of  10"1J  m/V) . 

The  signs  of  reXpti.  have  not  been  determined  except  where  specified. 


GaAs 

GaP 

Zinc  blende 
ZnSe  ZnS 

ZnTe 

CuCl 

ZnS 

Wurtzite 

CdS 

CdSe 

/  /  * 
c<c 

13.2 

12.0 

9.1 

8.3 

10.1 

7.5 

U 

C'* 

00 

9.5 

10.2 

/ 

-0.09 

-0.11 

-0.16 

-0.18 

-0.12 

-0.21 

-0.18 

-0.16 

-0.15 

tc*/e 

0.20 

0.23 

0.33 

0.35 

0.26 

0.27 

0.35 

0.41 

0.36 

+  1.03 

+  1.53 

*  2 .64 

+  2.93 

+  2.07 

-5.56 

+  3.63 

+  3.80 

+  3.61 

»Tvc 

-2.73 

-3.20 

-4.68 

-4.77 

-  6.41 

+  2.66 

-  5.63 

-6.71 

-7.40 

v  inco 

Mum 

-  1.7 

-  1.7 

-2.0 

-1.8 

-4.3 

-2.9 

-  2.0 

-2.9 

-3.8 

* 

-  1.5 

-  1.1 

2.0 

1.6 

4.3 

-2.4 

1.8 

2.4 

4.3 

*I«  P.  Kamlnow  and  E.  H.  Turner,  In  Handbook  of  Lasers,  edited  by  R.  J.  Pres¬ 
sley  (Chemical  Rubber  Co,,  Cleveland,  1971),  p.  453. 

bS.  Singh,  In  Handbook  of  Lasers,  edited  by  R.  J.  Pressley  (Chemical  Rubber 
Co.,  Cleveland,  1971),  p.  4H9. 

CI.  B.  Kobyakov,  Krlstallograflya  11,  419  (1966)  ISovlet  Phys.  Cryst.,  11,  369 
(1966)1. 


For  wurtzite  crystals,  we  neglect  the  small  dis¬ 
tortion  from  the  perfect  tetragonal  structure.  <9 

■  n 

can  be  expressed  in  terms  of  the  measured  macro¬ 
scopic  susceptibility  x  as  in  (7),  and  the  electro- 
optic  coefficients  of  zinc  blende  and  wurtzite 
crystals  are  obtained  as  follows:  zinc  blende, 

jrM1OD=0.3689a0:V/(ccVe);  (16) 

wurtzite, 

»'«,on=-2r13,on=0.4260aef(V/(^Ve),  (17) 

where  atfts=/3a02c0,  a0  and  c0  are  the  lattice  con¬ 
stants,  w  =  (e  -  1)(c<c-€)A2,  and  ru  are  in  units 
of  10' 13  m/V,  a0  and  aeff  in  angstroms.  Values  of 
the  parameters  edc' ,  /,  ec*  are  listed  in  Table  I. 

The  dependence  of  the  electro-optic  coefficients 
rm,on  on  the  bond  geometry  is  perhaps  the  most 
illuminating  feature  to  emerge  from  this  work. 


This  dependence  is  contained  in  curly  brackets 
in  (15).  For  diatomic  single-bond  crystals  is 
a  constant  and  the  geometrical  factor  becomes 

bonds 

D 

The  factor  /  is  typically  I/I  «  0.3.  Table  II  con¬ 
tains  a  listing  of  these  factors  for  some  key  di¬ 
rections  ( ijk )  in  crystals  of  the  zinc  blende,  wurt¬ 
zite  and  LiNbO,  classes.  It  follows  immediately 
that  when  Z^ni^C)  the  second  term  in  G  (Jk  is  an 
order  of  magnitude  larger  than  the  first  one.  In 
such  crystals  the  ionic  contribution  to  rilk  is 
about  an  order  of  magnitude  larger  than  the  elec¬ 
tronic  term.  Tliis  is  the  case  in  LiNb03  or  LiTaOr 
When  Zz„anj  =  0,  as  in  zinc  blende  and  wurtzite, 
we  have  to  settle  for  the  smaller  term  Z/^/,:i'n<aniQ;n»- 
This  is  the  main  reason  why  LiTa03  has  r33=  30.3 


TABLE  II.  Comparison  of  the  geometrical  factors  between  zinc  blende, 
wurtzite,  and  LINbO,.  The  quantities  shown  have  been  divided  by  the  num¬ 
ber  of  bonds  per  unit  cell,  nb  , 


Zinc  blende 

Wurtzite 

LINbOj 

(Nb-0),t,orI  (Nb-0)tonn 

0 

0 

0.475 

-0.669 

i 

4 

4 

4 

0 

*  t 

t 

0.184 

-0.185 

nrlZam  3 

0 

0.107 

-0.300 
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xl0'ls  m/V,  while  in  GaAs  and  ZnS,  r4l~2xlO"lJ 
m/V. 

The  effective  charge  e c*  is  related  to  the  Szigeti 
effective  charge  e,*  by  ec*  =  [(e„+ 2)/3e«,]e4*.  The 
value  of  e+/e  varies  from  0.2  to  0.4  and  seems 
to  be  independent  of  the  number  of  valence  elec¬ 
trons.  The  calculated  values  of  rUAlon  using  (17) 
and  (18)  are  shown  in  Table  I  as  rlonlc.  The  pure 
electronic  contribution  is  entered  as  reIec.  It  is 
obtained  directly  from  the  experimentally  deter¬ 
mined  second-harmonic-generation  coefficient  by 
rlk  =  -  4dlk/e2.  For  most  of  the  cyrstals  in  Table 
I,  yionic  is  positive  and  relec  is  negative.  There¬ 
fore,  the  predicted  electro-optic  coefficients, 
»'.um,hto=»'ionic+^eiec.  involve  the  algebraic  addi¬ 
tion  or  cancellation  of  two  numbers  of  compara¬ 
ble  values.  The  only  exception  is  CuCl.  Because 
of  the  unfilled  shell  in  Cu,  the  sense  of  bond  po¬ 
larization  in  CuCl  is  different  from  that  in  other 
crystals.7  The  signs  of  rlonlc  and  relec  of  CuCl 
are  thus  different  from  others.  However,  the 
magnitude  of  rlonic  is  larger  than  that  of  relec  in 
CuCl.  We  still  obtain  a  negative  electro-optic  co¬ 
efficient  for  CuCl.  The  predicted  values  in  Table 
I  are  In  good  agreement  with  experiment.  The 
worst  case  is  GaP.  It  is  interesting  to  note  that 
the  electronic  contribution  is  about  double  the 
ionic  contribution.  This  is  in  excellent  agree¬ 
ment  with  the  experimental  observation.8 

We  intend  to  extend  this  model  to  complex  crys¬ 
tals  with  different  point-group  symmetries.  The 
generalization  of  the  bond  parameters  used  above 


to  multibond  crystals  has  already  been  considered. 
The  one  parameter  which  will  need  added  thought 
is  e the  effective  ionic  charge.  It  was  found 
empirically  to  be  equal  to  C/Hw,  in  diatomic 
crystals.6  If  this  relation  a.ndec*  =[(e«,  +  2)/ 
3e„]e,*  are  valid  in  the  more  complex  crystals, 
then  our  model  can  be  applied  to  these  cases. 
Calculations  now  in  progress  on  KH2P04,  LiNb03, 
LiTaOj,  and  ternary  chalcopyrite  crystals  will 
be  reported  separately. 
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